from scipy.integrate import quad
import numpy as np

# 定义内部积分的被积函数
def inner_integral(t, x, n):
    return np.exp(-x * t) / t**n

# 定义外部积分的被积函数
def outer_integral(x, n):
    return quad(inner_integral, 1, np.inf, args=(x, n))[0]

# 计算积分
n = 2  # 你可以根据需要改变n的值
result, error = quad(outer_integral, 0, np.inf, args=(n,))
print(f"The result of the integral I_{n} is: {result}")